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We report imbibition experiments investigating the effect of evaporation 



f~^ ■ on the interface roughness and mean interface height. We observe a new ex- 

0^ . ponent characterizing the scahng of the saturated surface width. Further, we 



argue that evaporation can be usefuhy modeled by introducing a gradient in 



B 

i-rt ■ the strength of the disorder, in analogy with the gradient percolation model 

O ■ of Sapoval et al.. By incorporating this gradient we predict a new critical ex- 



ponent and a novel scaling relation for the interface width. Both the exponent 

value and the form of the scaling agree with the experimental results. 
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Recently the growth of rough interfaces has witnessed a veritable explosion of theoretical 
and experimental results, partly fueled by the broad interdisciplinary aspects of the subject 
||1|]. Much attention has focused on measuring the roughness exponent a, defined by the 
power law dependence upon the observation length scale I of the width w{t). Simulations on 
discrete models provide exponents in agreement with the predictions of phenomenological 
continuum theories [§. However, experimental studies find exponents significantly larger 
than the predictions of theory — for example, for dimension d = (1 + 1), theory predicts a = 
1/2 but experiments show a ~ 0.63 — 0.8 P,@|. Moreover, experimental studies frequently 
detect a crossover in w{i) to a different behavior above some characteristic length ix- It is 
currently believed that the anomalously large values of the exponent a is due to quenched 
pinning disorder |jl|J^,^ , however, a completely satisfactory explanation of the experimentally 
determined crossover length has not yet been found. 

Here we present imbibition experiments that probe the effect on the growth process of 
the evaporation rate and suspension concentration. We find that the scaling of the interface 
width changes with the evaporation rate and is characterized by a new exponent 7. We also 
present a model, inspired by the model of Ref. g], that predicts the experimentally-observed 
value of the new exponent characterizing the crossover effect. 

The key ingredient in the model is to allow for a gradient in the density of the pinning 
cells, which results in the stopping of the interface. Moreover, the detailed investigation 
of the scaling properties of the model provides us with a new scaling law and new critical 
exponents. Using this scaling law, we find good data collapse and scaling exponents that 
agree with the values determined analytically and numerically. 

In our experiments, paper — the "disordered medium" — is dipped into a reservoir 
filled with a colored suspension (coffee, ink) and the propagating wetting front is observed. 
The wetting front reaches a critical height, h^, above the level of the liquid, and stops 
propagating when the evaporation of the liquid induces the pinning of the interface by 
the inhomogeneities of the paper. We digitize the rough boundary between colored and 
uncolored areas and measure a roughness exponent a ~ 0.63 @]. 



Although the experiments are straightforward, their explanation in physical terms is less 
so. At microscopic length scales, paper is an extremely disordered substance, formed by long 
fibers that are randomly distributed and have random connections among one another. The 
wetting fluid propagates in these fibers mainly due to capillary forces, but the random nature 
of the fiber network and the particles in the suspension provide constant obstacles for the 
fluid flow 0. As we depart from the water source, evaporation is constantly decreasing the 
fluid pressure, making it more and more difficult for the fluid to overcome these microscopic 
"obstacles". At the critical height, the fluid pressure balances the effect of the pinning 
obstacles and the fluid stops propagating 0]. 

We anticipate, on physical grounds, that the smaller the evaporation, the larger the 
critical height will be. To check this intuition we repeated our experiments in environments 
with different rates of evaporation. As we decreased the evaporation rate, the height reached 
by the interface increased sharply. These effects can be observed in the pinned interfaces 
obtained experimentally (Fig. |l|(a)). 

While the above results can be understood from microscopic considerations, the effect 
of the evaporation or of the suspension density on the scaling of the interface roughness is 
nontrivial. 

To understand these results, we propose a model to describe the formation of the pinned 
interface in a disordered medium. In (1 + 1) dimension, we model the pinning obstacles by 
blocking, in a lattice of horizontal size L, a fraction p{h) of the cells in each horizontal row, 
where h is the height from the bottom of the lattice. We start, at t = 0, from a horizontal 
line of wet cells at the bottom edge of the lattice. At time t + 1 we wet all unblocked cells 
which are nearest neighbors to the wet region at time t. We also apply the rule that every 
cell, blocked or not, bellow a new wet cell becomes wet as well (Fig. ^). The motivation 
for this rule is the experimental observation that the wet region is, at least at macroscopic 
length scales, nearly free of dry islands. 

If p{h) = Po, this model generates an interface which propagates with a constant {p- 
dependent) velocity ii Po < Pc — 0.47, and becomes pinned by a directed percolation cluster 



that spans the system at po > Pc @jl]- However, although the actual disorder in the 
paper is not height dependent, its effect in pinning the propagation of the fluid is increasing 
with height, due to the decrease in the fluid pressure. The most physical assumption is an 
exponential decrease of the fluid pressure or, equivalently, of the driving force. This will 
lead to an "effective" increase in the density of pinning obstacles as we depart from the 
reservoir, i.e., p = p{h). Hence 

p(/i)-PoOcl-e-^/'^°. (1) 

If h <^ ho and Pc — Po ^ Pc, we can write 

p{h) - Po oc h-^h oc {Vp)h. (2) 

Hence, in this limit, we find a constant non-zero gradient in the density of pinning obstacles. 
The presence of the gradient Vp changes the width of the pinned interface (Fig. 0(b)) and 
its scaling form (Fig. ^). Our simulations show that for observation scales i much smaller 
than some characteristic crossover length i^, the saturated width behaves as w ^ £", but 
for i':^ ix, the width saturates at a value Wsat that depends upon the gradient as 

Ws.t ~ (Vp)'^. (3) 

This behavior can be expressed by a scaling law of the form 

«;(£,Vp)~r/(£/£x) (4a) 



,^ r^ I x/nl-'r/" 



(Vp)-^/". (4b) 



The scaling function f{u) satisfies f{u^ 1) ~ const and /(m ^ 1) ~ m °'. Our simulations 
(see Fig. ^ for a system of size L = 16384 yield the exponents 

asim = 0.63 ± 0.02, 7sim = 0.52 ± 0.02. (5) 

We remark that the validity of the scaling law (^ and the values of the exponents do not 
depend on the exact form oi p{h) but only on the value of Vp{h) at h^ [§]. 



The value of a can be understood from the mapping to directed percolation, since the 
conditions for a complete pinning of the interface do not change from the models of Refs. [^Q . 
In directed percolation, the size of a cluster is characterized by a longitudinal correlation 
length ^11 and a transverse correlation length ^_l that, near pc, behave as 

^11 ~ be - pI"''" , i^-^\vc-v\~''^- (6) 

The roughness exponent is related to the exponents of directed percolation as 0,|] 

a = z/j_/z/||. (7) 

Using the known values of v^ and v\\_ in relation ^ we predict a = 0.633 ± 0.001, in 
agreement with our simulation result (^. 

The exponent 7 can be related to ux theoretically. A point of the interface, at distance 
Wsat of the critical height, is pinned by a directed percolation cluster if the transverse size 
of that cluster is of order ^±{p). At that point we have p = p{hc ± Wsat) ~ Pc ± u^sat ■ Vp. 
Therefore, using Eq. (^ we find |8|JT0[ 



Wsat ~ ^±{P) ~ \Pc - {Pc ± M^sat " Vp) | ''^ , 

U-sat-ksaf Vp]-'^^. (8) 

From Eqs. (|) and (|^) follows 

7 = z/x/(l + z/±). (9) 

Since i/_l is known accurately 0, Eq. (^ predicts 7 = 0.523 ± 0.001, in excellent agreement 
with our simulation result (|). 

Without the gradient, the interface has critical behavior only if we tune p to pc- However, 
with the gradient the interface always stops at the critical height h^. This critical height 
can be calculated from the condition p{h^ = p^. Thus from (^ we obtain 

he ~ {Vp)-\ (10) 



i.e. the height reached by the wetting fluid is inversely proportional to the gradient in the 
disorder. 

The experimental data presented in Fig. ^(a) remarkably resemble the data obtained for 
the model. However, without knowing the actual value of the gradient in the experiments, 
it is not possible to check the validity of the scaling law (^ experimentally. Nonetheless, 
measuring the critical height in the experiments and using Eq. ([10|) , we are able to estimate 
Vp, the gradient in the "effective disorder" , for the experiments, up to a multiplicative con- 
stant. Using these experimentally determined values of Vp, we rescale the results obtained 
for the width according to the scaling law @. In Fig. ^(b) we show this rescaling, where 
we used 

Oexp = 0.65 ± 0.05, 7exp = 0.49 ± 0.05, (11) 

The experimental values of both exponents agree well with the results obtained from the sim- 
ulations (Table |I|) and with the theoretical predictions based in known results from directed 
percolation. 

The generalization of the model to d = (2 + 1) is straightforward, and in this case the 
pinned interface can be mapped to directed surfaces |^ , a percolation problem that has not 
been thoroughly investigated. We simulated the model for a 512 x 512 system; the critical 
exponents that give the best data collapse are 

ttsim = 0.43 ± 0.04, 7sim = 0.32 ± 0.02. (12) 

From these results, we calculate the exponents characterizing the transverse and longitudinal 
correlation lengths for the directed surfaces problem, obtaining 

u± = 0A7±0M, i^ij = l.liO.l. (13) 

In summary, we have performed imbibition experiments to study the effect of evaporation 
on interfacial phenomena. We have also developed a model that incorporates evaporation 
by introducing a gradient in the density of pinning cells [0 . The model provides insight into 
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three previously unexplained aspects of imbibition experiments: (i) The interface always 
stops growing, after some finite time. Due to the gradient, the wetting interface only moves 
until it reaches a critical concentration of pinning cells. This gradient in pinning cells arises 
from the balance between the evaporation of the fluid and the capillary forces tending to 
move it along the paper, (ii) The flnal height of the interface, he, increases when the 
evaporation is reduced, due to the smaller effective gradient in the pinning disorder, (in) A 
new exponent 7 was found characterizing the dependence on the gradient of the saturation 
width and the characteristic length i^- Good agreement was found between experimental, 
analytical and simulation values of the exponents. 

This work is part of the Ph. D. thesis of L. A. N. Amaral and A.-L. Barabasi. We 
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National Science Foundation. 
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FIGURES 
FIG. 1. Photographs of pinned interfaces in: (a) Imbibition experiments with coffee and 

paper towels for (i) high evaporation rate: {'Vp)exp = 0.94(/o, and (ii) low evaporation rate: 

(Vp)exp = 0.25go- Here go is the undetermined multiplicative constant discussed in the text. 

(b) Simulations of the model, with L = 256, for different values of the gradient: (i) Vp = 2^^ and 

(ii) Vp = 2^^^. Readily apparent from these photographs is the increase in both the final heights 

and widths of the interface with the decrease of the gradient. 

FIG. 2. Example of the time evolution of the model for a very small lattice (L = 5). Here, grey 
squares represent blocked cells and white squares represent unblocked cells. The numbered cells 
are wet. The numbers indicate at which time step the cells first become wet. At t = 4, we wet the 
cells at the left and at the right of the cell numbered 3. Also, in the same time step we wet the 
cells below those two, regardless of the fact that they were previously blocked. Similarly, at t = 5, 
we are able to wet cells in the first column from wet cells in the second column that were, at some 
earlier time, blocked cells. The heavy line indicates the pinned interface. 

FIG. 3. The simulation results for the width w{i, Vp) of the pinned interface, (a) The widths 
for several values of the gradient (averaged over 512 runs for each value of the gradient), (b) The 
same simulation results, plotted in the scaling form of Eq. (3), using the values of the exponents 
from Eq. (4). 

FIG. 4. The experimental results for the width w{i, Vp) of the pinned interface, (a) The widths 
for several values of the gradient (in units of go). The values of the gradients were calculated as 
described in the text; the error in these values is smaller than 10%. The widths were corrected 
by a multiplicative factor to make them coincide for the smallest i. (b) The same experimental 
results, plotted in the scaling form of Eq. (3), using the values of the exponents from Eq. (10). 
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TABLES 
TABLE L Critical probability and exponents for dimension (1 + 1) and (2 + 1), calculated from 

the simulations of the present model. 

dim. (1 + 1) (2 + 1) 

Pc 0.47 ±0.03 0.75 ±0.03 

a 0.63 ±0.02 0.43 ± 0.04 

7 0.52 ±0.02 0.32 ±0.02 

u_i 1.09 ± 0.08 0.47 ± 0.04 

Ml 1.7 ±0.1 1.1 ±0.1 
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